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We calculate the neutrino mean free path in color superconducting quark matter, and employ it to 
study the cooling of matter via neutrino diffusion in the superconducting phase as compared to a free 
quark phase. The cooling process slows when quark matter undergoes a second order phase transition 
to a superconducting phase at the critical temperature T c . Cooling subsequently accelerates as the 
temperature decreases below T c . This will directly impact the early evolution of a newly born 
neutron star, should its core contain quark matter. Consequently, there may be observable changes 
in the early neutrino emission which would provide evidence for superconductivity in hot and dense 
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In this article we study heat diffusion via neutrinos in dense, color superconducting quark matter. Recent theoretical 
works jlp^l suggest that quarks form Cooper pairs in medium, a natural consequence of attractive interactions 
destabilizing the Fermi surface. This would likely affect the early evolution of neutron stars born through a Type 
II supernova explosion, where the central role of neutrino diffusion through a strongly-interacting medium has long 
' been postulated on theoretical grounds ||. 

, Type II (core collapse) supernovae are triggered by the implosion of the inner core of a massive star (M sta r ^8 — 20 

ly-j ■ M ), when the core mass is on the order of the Chandrashekar mass (M core ~ 1.4 M Q ). During the implosion nearly 
all (~ 99%) of the enormous gravitational binding energy (~ 10 53 ergs) gained is stored as internal energy of the newly 

C^) ■ born, proto- neutron star (PNS). The subsequent evolution of the proto-neutron star is driven by neutrino diffusion. 
Temporal and spectral characteristics of the neutrino emission depend on the rate at which neutrinos diffuse through 
the imploded PNS which, at this early stage, is composed of hot (T ~ 20 — 30 MeV) and dense (ns ~ 2 — 3 no where 
O |' no — 0.16 fm -3 ) strongly-interacting matter. 

Neutrino emission during this phase is a directly observable feature of a galactic supernova explosion. The few 
neutrinos detected from SN 1987A indicate that neutrinos of mean energy (E„) ~ 20 MeV are emitted on a diffusion 
time scale of about 10 — 20 seconds. It is reasonable to expect that neutrino mean free path in the inner, denser 
regions of the star will strongly influence the temporal characteristics of a supernova neutrino signal. Supernova 

. £h ' neutrinos can therefore reveal properties of matter at high baryon density, at temperatures substantially lower than 

^ . those expected in relativistic heavy ion collisions. 

Although the idea of quark pairing in dense matter is not a new one 0] , it has recently seen renewed interest in 
the context of the phase diagram of QCD Q. Model calculations, mostly based on four-quark effective interactions, 
predict the restoration of spontaneously broken chiral symmetry through the onset of color superconductivity at low 
temperatures ||] . For much higher densities color superconductivity is manifest through perturbative gluon exchange 
HQ, which can be calculated systematically suggests that the phenomenon is robust. For densities and temperatures 
relevant to neutron stars, quark matter is therefore expected to be superconducting. Models generally predict an 
energy gap of A ~ 100 MeV for a typical quark chemical potential of a q ~ 400 MeV. As in BCS theory, the gap 
will weaken for T > and at some T = T c there is a (second-order) transition to a "standard" quark-gluon plasma. 
Thus, when cooling from temperatures greater than critical, the formation of a such a gap in the fermionic excitation 
spectrum in quark matter at high density will influence various neutron star observables - if neutron stars contain 
quark matter in their cores at early timdj. Two examples recently investigated are the effects on magnetic fields Jicf] 
and on the thermal evolution of neutron stars at late time pT| ], when interior temperatures evolve from T ~ 1 MeV 
to T< 1 KeV. In this work we shall consider the cooling of quark matter at earlier times, when temperatures pass 



This remains unclear. In particular, the existence of a finite electron neutrino chemical potential at very early time has 
shown to inhibit the appearance of quark matter. See Ref. M for a review. 
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through the critical T c ~ A. Assuming a simplified scenario, we investigate the thermal evolution of the inner core of 
a proto-neutron star as it cools via neutrino diffusion during its first few seconds. 

Our main finding is that the neutrino mean free path in the superconducting phase has a strong dependence on 
temperature. An energy gap tends to increase the neutrino mean free paths exponentially when T -C A. However 
in the intermediate regime, when T ~ A, the temperature dependence is not exponential. From this we predict a 
uniquely uneven cooling process for simplified PNS matter, marked by a slowdown of the cooling rate when the system 
undergoes a second-order phase transition. This is a consequence of the specific heat being peaked at T c , a standard 
characteristic of a phase transition, rather than a large change in the neutrino mean free path. Since the energy gap 
vanishes at this point, the mean free path is not significantly modified in the neighborhood of T c . 

In Section II the neutrino mean free path, always denoted A in this work, is computed in a background of supercon- 
ducting quarks. Taking a general, BCS-type model for the energy gap, we relate A to in-medium quark polarizations 
via the differential and total neutrino-quark cross sections. Explicit formulae are derived for the imaginary parts of 
the vector and axial-vector polarizations of a relativistic Fermi system with an energy gap. After assuming BCS-type 
mean field behavior of the gap and the specific heat, in Section III we compute a characteristic time for heat diffusion 
via neutrino emission from a simple model of superconducting quark matter. We then consider some astrophysical 
consequences and outline the model's applicability to proto-neutron stars. Section IV contains our conclusions. 



II. NEUTRINO-QUARK SCATTERING IN A COLOR SUPERCONDUCTOR 

The primary process by which heat escapes a proto-neutron star is neutrino diffusion, and so a significant con- 
sequence of color superconductivity in this context will be modified neutrino propagation. While noting that the 
neutrino production rate will also differ from that of normal matter, in the diffusive regime one can see that the 
dominant critical behavior will be a change in the inelastic quark-neutrino cross section, since here neutrino produc- 
tion rates decouple from the transport equation and depend only on the neutrino mean free path. In this section we 
calculate the differential and total cross sections and then compute the neutrino mean free path in two-flavor quark 
matter. The magnitude of the superconducting gap, A, is taken as arbitrary within a range of values found in recent 
literature. Closely following BCS theory, we assume the gap to be a constant, and calculate the response functions 
and neutrino cross sections in the weak coupling approximation. 

The neutral current coupling between neutrinos and quarks, a four-fermion effective interaction for energies E v <C 
Mz, is written as 

G F 

£w = -j= ^7m(! ~ l5)vq^{C v - C A 75)q, (1) 

where Gf = 1.166 x 10~ 5 GeV~ 2 is the Fermi weak coupling constant and Cy and Ca are the flavor-specific vector 
and axial vector coupling constants, respectively. The differential neutrino scattering cross section per unit volume in 
an infinite and homogeneous system of fermions as calculated in linear response theory is [T^JT^ | 

1 #o g| E 3 [l-f„(E 3 )] ( a0 , m 

V<PSl 3 dE 3 ~ 32-n- 2 Ei [1 - exp (—q /T)] Q/3j ' [ ) 

where Ei (E 3 ) is the incoming (outgoing) neutrino energy. The factor [1 — exp(— qo/T)]^ 1 maintains detailed balance 
and the final state blocking of the outgoing neutrino is enforced by the Pauli blocking factor, [1 — f v {E 3 )\. The 
neutrino tensor L a p is given by 

L af3 = %[2k a k P + (k ■ q)g al3 - (fcV + q 01 ^) T ie a ^ u k^q u ] , (3) 

where the incoming four-momentum is k a and the momentum transferred to the medium is q a . The minus (plus) 
sign on the final term applies to neutrino (anti-neutrino) scattering. 

The medium is characterized by the quark polarization tensor IIq,^. In the case of free quarks, each flavor contributes 
a term of the form 

IM<?) = -*Tr c J Tr [So(p)r a S (p + q)Tp] , (4) 

where So(p) is the free quark propagator at finite chemical potential and temperature. The outer trace is over color 
and simplifies to a N c = 3 degeneracy. The inner trace is over spin, and the T a are the neutrino-quark vertex functions 
which determine the spin channel. Specifically, the vector polarization is computed by choosing (r^r^) = (70,7,3). 
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The axial and mixed vector-axial polarizations are similarly obtained from (ra,^) = (7075,7/375) and (r„,r^) = 
(7a:7/375), respectively. 

The free quark propagators in Eq. (^]) are naturally modified in a superconducting medium. As first pointed out 
by Bardeen, Cooper, and Schricffer several decades ago, the quasi-particle dispersion relation is modified due to the 
presence of a gap in the excitation spectrum. In calculating these effects, we will consider the simplified case of QCD 
with two quark flavors which obey SU(2)^x SU(2)# flavor symmetry, given that the light u and d quarks dominate 
low-energy phenomena. Furthermore we will assume that, through some unspecified effective interactions, quarks pair 
in a manner analogous to the BCS mechanism jl4| . The relevant consequences of this are the restoration of chiral 
symmetry (hence all quarks are approximately massless) and the existence of an energy gap at zero temperature, Ao, 
with approximate temperature dependence, 



A(T)=A yi-^-j . (5) 

The critical temperature T c ~ 0.57Ao is likewise taken from BCS theory; this relation has been shown to hold for 
perturbative QCD Jig ] and is thus a reasonable assumption for non-perturbative physics. 

Breaking the SU C (3) color group leads to complications not found in electrodynamics. In QCD the superconducting 
gap is equivalent to a diquark condensate, which can at most involve two of the three fundamental quark colors. The 
condensate must therefore be colored. Since the scalar diquark (in the 3 color representation) appears to always be 
the most attractive channel, we consider the anomalous (or Gorkov) propagator 

F(p)abfg = (qj a (p)Cj 5 q gb {-p)) 



= -^e/ ff A(^| + ^|]7 5 C. (6) 



Here, a, b are color indices, /, g are flavor indices, e a bc is the usual anti-symmetric tensor and we have conventionally 
chosen 3 to be the condensate color. This propagator is also antisymmetric in flavor and spin, with C = —17072 being 
the charge conjugation operator. 

The color bias of the condensate forces a splitting of the normal quark propagator into colors transverse and parallel 
to the diquark. Quarks of color 3, parallel to the condensate in color space, will be unaffected and propagate freely, 
with 

So( P t f = i* a S} + (p M y. _ M70 ) . (7) 

Y-Fo -^p Po p J 

This is written in terms of the particle and anti-particle projection operators A + (p) and A~(p) respectively, where 
A ± (p) — (1 ± 707 • p)/2. The excitation energies are simply E p = \p\ — \i for quarks and E p = \p\ + /j, for anti-quarks. 

On the other hand, transverse quark colors 1 and 2 participate in the diquark and thus their quasi-particle propa- 
gators are given as 



*<ws = m (£4 + pi ) (p,f - ™> . (8) 

¥0 Sp Po Sp 



The quasi-particle energy is £ p — (|p| — /1) 2 + A 2 , and for the anti-particle £ p — (\p\ + p) 2 + A 2 . 

The appearance of an anomalous propagator in the superconducting phase indicates that the polarization tensor 
gets contributions from both the normal quasi-particle propagators (j^) and anomalous propagator (^). Thus, to order 
G 2 ?, Eq. §Q) is replaced with the two contributions corresponding to the diagrams shown in Fig. [l], and written 

n Q /3(q) =-ij {Tr [So(p)T a S (p + q)Tp] + 2Tr [S(p)T a S(p + q)T p ] + 2Tr [F(p)T a F(p + q)Tp]} . (9) 

The remaining trace is over spin, as the color trace has been performed. Fig. ^(a) corresponds to the first two terms, 
which have been decomposed into one term with ungapped propagators and the other with gapped quasi-particle 
propagators (||). Fig. |l|(b) represents the third, anomalous term. 

For neutrino scattering we must consider vector, axial, and mixed vector-axial channels, all summed over flavors. 
The full polarization, to be used in evaluating Eq. (||), may be written 
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(a) (b) 

FIG. 1. Standard loop (a) and anomalous loop (b) diagrams contributing to the quark polarization operator. 

u a0 = J2 [{c f v ?K P + (c f A fu^ - 2C v c{n v j] . (io) 
/ 

The coupling constants for up quarks are C v = \ — | sin 2 9w and G\ — \ , and for down quarks, C v = — \ + 1 sin 2 9w 
and C d A = —\, where sin 2 6w — 0.23 is the Weinberg angle. 

If we specify a frame in which the transfer momentum is — (go, 9, 0, 0) we can separate longitudinal components 

as 



LL L — ~2 n oo — ~2 LL n — ~ U 10 ' \ LL ) 



- 2 a 2 a 2 

~ 2 11 — 

% qoq 

and the transverse, 

n£ = n 2 ^ 2 = n 3 v 3 . (12) 

Identical definitions apply to the axial polarizations and n^. The non-zero mixed correlation function is written 

n v J = ie a ^ q »n VA . (13) 

Detailed calculations of these polarization functions are given in the Appendix. 

As the gap increases, the superconducting quasi-particles naturally become the dominant excitations of the back- 
ground, a property clearly visible in the neutrino response functions. The left panel of Fig. ^| shows the longitudinal 
response in the vector channel. The free quark case, shown as a solid line labeled A — 0, is the standard result describ- 
ing Pauli-blocking and kinematics of massless single particle excitations [ jOf . Here, energy- momentum conservation 
restricts the response to the spacelike region (q < q). Superconductivity modifies this result, as the quasi-quark 
excitations become suppressed due to the pairing correlations at the Fermi surface. At the same time, the response 
is enhanced when qo > 2A, signifying the excitation of Cooper pairs. In particular, the results with greater A clearly 
show the threshold for these excitations at energy transfer qo — 2A. Results for A = 10, 30 and 50 MeV show the 
gradual reallocation of response strength from small qo to the region qo > 2A. Since scattering probes only the 
spacelike region, the contribution to the cross section is generically suppressed in the superconducting phase. 
Contributions at q^ > q will contribute to the neutrino production rate, rather than scattering cross section, when 
the temperature is near A. Analytic results may be obtained for small transfer energies. When go is much smaller 



than all other energy scales, the vector longitudinal response (see Eq. (Al) of the Appendix) reduces to 



Km Im Il£(A) = — 1— n^(A = 0) 

go^o 1 + e pzi 

,2 



2nq 1 + eP* ' 1 ' 

From this the weakening of the low energy, vector-longitudinal response can be calculated for a given gap A. 

The axial-longitudinal response, which physically corresponds to the excitation of spin waves, is shown in the right 
panel of Fig. S. As with the vector channel a threshold of 2 A is apparent but, unlike the previous case, the response 
as qo — > is enhanced. This is manifest in the the limit qo <C T ~ A, where one finds jl7| 

Im Uf (A) ~ A sech 2 In (^j n£ (A = 0) , (15) 
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FIG. 2. Vector-longitudinal (left) and axial-longitudinal (right) responses as a function of energy transfer. The momentum 
transfer is fixed at q = 50 MeV. 



where in this limit n^(A = 0) = [i 2 / Airq. This logarithmic enhancement will lead to an integrable peak in the 
differential cross section at qo = 0, as will be described below. 

The transverse response functions, 11^ and 11^ as defined in Eq. (|l^), exhibit behavior similar to the longitudinal 
channels. The primary distinction is that superconductive coherence at low qo enhances the vector channel and 
suppresses the axial. In the former case the interference effects are constructive and in the latter case they are 
destructive |Q, which is a reversal of the longitudinal results. This difference has been theoretically understood and 
experimentally observed in electric superconductors, where the absorption of acoustic waves (vector-longitudinal) is 
suppressed while the infra-red absorption (vector-transverse) is enhanced for small energy transfer p9| , p0| p[ 

In addition to the vector and axial polarizations, the mixed vector-axial polarization (|l^) also contributes to the 
cross section. While this contribution is always much smaller - by at least one order of magnitude - we note that it is 
enhanced in the superconducting phase, leading to an amplified difference between neutrino and anti-neutrino cross 
sections in PNS matter. 

Once the polarization tensors have been computed it is straightforward to obtain the differential cross section, 
Eq. (||) . Results for neutrinos of energy E v = 50 MeV, ambient matter conditions of fj, q = 400 MeV and T — 30 MeV, 
and gaps of varied size are plotted as a function of transfer energy in Fig. ^. A striking feature of these results is 
the singular behavior of the differential cross section near go = in superconducting matter. This is the (integrable) 
logarithmic divergence of Il£ and (see (Eq. [H])). The threshold behavior seen in differential cross section, for the 
case A = 10 and 30 MeV, at qo = 2A and qo = — 2A correspond to excitations of the Cooper pairs. 

The total cross section (per unit volume) is obtained by integrating over all neutrino energy transfers and angles. 
From this the mean free path is determined, since 

» -(£)"'■ < M » 

As suggested by the differential cross section in Fig. g, the total cross section is reduced in the presence of a gap A. 
The logarithmic peak at qo = has a minimal effect after integration, when A ~ T and when the neutrino energy 
E v ~ ttT since neutrino probe a significant part of the response outside of this qo — region. 



2 There is naturally a vast literature concerning coherence effects in superconducting metals, and we find it both gratifying 
and reassuring to obtain similar results in this less terrestrial context. 
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FIG. 3. Differential cross section as a function of energy transfer and for neutrino energy E v — 50 MeV. 



Results for the neutrino mean free path, A, are shown in Fig. [l| as a function of incoming neutrino energy E v 
(ambient matter conditions of [i q — 400 MeV and T — 30 MeV have again been used). The same energy dependence 
has been computed previously in the case of free relativistic and degenerate fermionic matter it decreases as 
l/E% for E v T, and 1/E V at E v <C T. The results indicate that this energy dependence is not modified by the 
presence of a gap when A ~ T . Thus the primary effect of the superconducting phase is a much larger mean free 
path. This is consistent with the suppression found in the vector-longitudinal response function, which dominates the 
sum polarization sum (|l0|), at q < q. 



III. THE COOLING OF AN IDEALIZED QUARK STAR 



The immediate application of the last section's results is neutrino emission from a proto-neutron star containing 
quark matter. While we will fall far short of a complete description of the role of color superconductivity in such 
a complicated environment, in this section we will outline the key ingredients for estimating the possible physical 
observables. The first subsection describes a simple and general model for the cooling of quark matter through 
temperatures relevant to proto-neutron stars, and the second addresses the applicability of this scenario to more 
realistic situations. 



A. The Cooling of Superconducting Quark Matter 

Having determined the neutrino mean free path in a color superconductor, we now consider the cooling of a 
macroscopic sphere of quark matter as it becomes superconducting. As stated previously, this toy model is motivated 
by the possibility that the core of neutron star might contain such matter. Following our preceding calculations, we 
will consider the relatively simple case of two massless flavors with identical chemical potentials. Furthermore, we 
will disregard the quarks parallel in color to the condensate; i.e. we consider a background comprised exclusively of 
quasi-quarks. 

The cooling of a spherical system of quark matter from T ~ T c ~ 50 MeV is driven by neutrino diffusion, for the 
neutrino mean free path is much smaller than the dimensions of system of astrophysical size, and yet several orders 
of magnitude larger than the mean free path of the quarks. The diffusion equation for energy transport by neutrinos 
in a spherical geometry is 

tfy— = --— , (17) 
dt r 2 dr ' 
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FIG. 4. Left Panel: Neutrino mean free path as a function of neutrino energy. Right Panel: The neutrino mean free path for 
neutrino energy E v = ttT, plotted as a function of the gap A. Results, shown here for T — 30 MeV, are virtually independent 
of temperature for T< 50 MeV. 



where Cy is the specific heat per unit volume of quark matter, T is the temperature, and r is the radius. The neutrino 
energy luminosity for each neutrino type, L u , depends on the neutrino mean free path and the spatial gradients in 
temperature and is approximated by 

L v = -6 J dE v E yX(E v )^^- , (18) 

where c denotes the speed of light in vacuum. In our analysis we assume that neutrino interactions are dominated 
by the neutral current scattering common to all neutrino types. Consequently, we take the same neutrino and anti- 
ncutrino mean free path for every neutrino flavor, giving rise to the factor of six in Eq. (|f). The equilibr ium Fermi 
distribution, f(E„), and the (scattering) mean free path, \{E U ) 1 are integrated over all neutrino energies, E v . 

The solution to the diffusion equation will depend on the size of the system and its initial temperature gradients. 
However, we are merely interested in a qualitative description of cooling through a second-order phase transition to 
superconducting matter. The temporal behavior is characterized by a time scale r c , which is proportional to the 
inverse cooling rate and can hence be deduced from Eq. (0). The characteristic time 

r c (T) = C v (T)^—, (19) 

is a strong function of the function of the ambient matter temperature since it depends on the matter specific heat and 
the neutrino mean free path. This applies to a system characterized by the radial length R and the energy-weighted 
average of the mean free path, (A(T)). Following our general treatment of the superconducting gap, we assume that 
the temperature dependence of the specific heat is described by BCS theory. We will then use the results obtained in 
the previous section to calculate (A(T)). Furthermore, since neutrinos are in thermal equilibrium for the temperatures 
of interest, we may assume 

(A(T)) ~ \{E V = ttT) . (20) 

The quantity on the left is dependent on the gap A, a dependence computed in the previous section and plotted 
in the right panel of Fig. ^. The results indicate that for small A/T the neutrino mean free path is not strongly 
modified, but as A/T increases so too does A, non-linearly at first and then exponentially for A/T> 5. We note that 
the diffusion approximation is only valid when A <ti R and will thus fail for very low temperatures, when A ~ R. 
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FIG. 5. The figure shows the extent to which different physical quantities are affected due to the superconducting transition. 
Ratios of the cooling time scale (solid curve), the inverse mean free path (short-dashed curve) and the matter specific heat 
(dot-dashed curve) in the superconducting phase to that in the normal phase is shown as a function of the matter temperature. 
The ratio of the gap to it zero temperature value Ao is also shown (long-dashed curve). The quark chemical potential is 
fi q = 400 MeV and A = 100 MeV. 

The ratio t*(T)/t c (T), a measure of the extent to which the cooling rate is changed by a gap, is shown by the solid 
line in Fig. |^. The ratio A/A A , plotted with the short-dashed curve, measures the decrease in neutrino interaction 
rates in the superconducting background. The temperature dependences we have taken from BCS theory, that of 
the specific heat (dashed curve) and the magnitude of the gap itself (dot-dashed curve), are shown for reference in 
Fig. ^|. The results shown in Fig. |^ are readily interpreted. The cooling rate around T c is influenced mainly by the 
peak in the specific heat associated with the second order phase transition, since the neutrino mean free path is not 
strongly affected when A«T. Subsequently, as the matter cools, both Cy and A -1 decrease in a non-linear fashion 
for A ~ T . Upon further cooling, when A ^> T, both Cy and A -1 decrease exponentially. Both of these effects 
accelerate the cooling process. 

We conclude that if it were possible to measure the neutrino luminosity from the hypothetical object described here, 
a unique temporal profile would be observed. This suggests that if a second order, superconducting transition were 
to occur inside a proto-neutron star it could be identified by the temporal characteristics of the late time supernova 
neutrino signal. Specifically, there would be a brief interval during which the cooling would slow around T ~ T Cl 
signified by a period of reduced neutrino detection. 

B. Relation to Proto-Neutron Stars 

The temporal pattern of neutrino emission deduced here would be observable evidence of the onset of color supercon- 
ductivity in dense matter. However, ours is a simple model, and we must temper our speculations with considerations 
of more realistic systems. 

Neutrinos emitted from the core of a proto-neutron star, be it hadronic or quark matter in any phase, must pass 
through a large amount of matter in the outer shell. The neutron-rich material is characterized in the outer shell is 
opaque to neutrinos and will thus even out any sharp temporal features associated with the interior emission. This 
is the first and foremost of caveats since this will directly impact the possibility of observing a characteristic in the 
neutrino signal associated with the phase transition. 

As discussed in Section II, the scalar diquark condensate in a two-flavor color superconductor is necessarily colored. 
Therefore quarks of one color, taken conventionally as color 3, will be color-orthogonal to the scalar condensate and 
can remain ungapped. We have not included these color-3 quarks in our analysis of neutrino scattering since their 
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fate is uncertain |2l]] : they could form color-symmetric Cooper pairs, or perhaps a chiral condensate. If we take the 
simplest scenario and assume that they remain free, their presence in the medium will further dilute any direct effects 
of the gapped quasi-quarks. Specifically, the physical quantities in plotted in Fig. [5] will not vanish when T — > 0, 
instead being reduced to one third of their A = values. Likewise, the maxima at T c will be reduced relative to a 
color-neutral background. 

The second-order phase transition is taken directly from BCS theory, which we assume rather than derive. While 
this is the mean- field result for two massless flavors with equivalent chemical potentials, neutron star matter is 
constrained in two ways. First of all, weak- interaction equilibrium requires fid — fx u = fi e , where the subscripts refer 
to down quarks, up quarks, and electrons. But since a finite electron number is required to achieve electric charge 
neutrality in the star, /i e cannot vanish and we necessarily have fid 7^ This difference in chemical potentials is 
likely to drive a first order rather than a BCS-second order transition p^ ]. 

The other notable omission in our study is the strange quark. As before, weak interaction equilibrium requires that 
fi s = fid, and thus for m s strange quarks will be present. Furthermore, the finite strange quark mass implies a 
mismatch in the Fermi momenta which would also drive a first order transition Q . A generic consequence of a first 
order transition is a mixed phase containing both normal and superconducting quark matter, and transport in the 
heterogeneous mixed phase is qualitatively different from that considered here, for neutrino scattering will depend on 
the size and nature of the structures (droplets) present. In previous work it was shown that neutrino mean free path 
in the heterogeneous phase can be greatly reduced due to coherent scattering of droplets in the mixed phase p3| . 
Combining the results of Section II with neutrino transport in a mixed phase of superconducting and normal quark 
matter is beyond the scope of this work. 

Therefore, while our calculation of the neutrino mean free path will be an essential ingredient in a more realistic 
and hence more complicated model of neutron star evolution, our toy model only applies to a highly idealized quark 
core of a neutron star where A> \fid — fi u \- 

IV. CONCLUSIONS 

Motivated by the physical relevance of neutrino diffusion in the cooling of strongly-interacting matter, we have 
analyzed the effects on neutrino-quark scattering arising from a second-order phase transition from normal to color 
superconducting quark matter. The principal microscopic ingredient is the neutrino mean free path, and this was 
calculated in linear-response theory with a BCS-type superconducting background. We then applied this result to a 
schematic model of quark matter at temperatures and densities relevant to the evolution of proto-neutron stars. The 
modified mean free path for neutrinos in a superconducting background was computed in the simplified approximation 
of iso-symmetric, two-flavor quark matter. We have enumerated the main shortcomings of these simplifications and 
realize that a more realistic treatment of PNS evolution would include many other effects of similar importance. 

Despite these caveats, our toy model calculation indicates that a superconducting transition in the quark core of a 
PNS can impact its early cooling and thereby potentially alter the temporal characteristics of the neutrino emission. 
We view this work as a first step towards an understanding of how the presence of color superconducting matter in 
the core of a neutron star may affect the early - and observable - neutrino signal. Quark pairing invariably occurs in 
theoretical treatments of finite-density QCD, and thus our microscopic calculations of the neutrino cross sections are 
pertinent to transport processes in dense quark matter. Given the real prospect of detecting neutrinos emitted from 
a future supernova event, such transport processes might someday serve as an probe of the properties of extremely 
dense matter. 
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APPENDIX: QUARK POLARIZATION TENSORS 



At densities relevant to this work, only the quark particle-hole excitations are accessible. Thus we discard all 
anti-particle and anti-hole contributions from the propagators in Eqs. (||) and (||). 
The imaginary part of the vector longitudinal polarization is, for each flavor, 



Im n^(g , q) = -2ir 



d 3 p d 3 k -o. , W1 „ f . 

■S 3 (p- k - q)(l + p- k) 



(2tt) 3 (2tt) 3 



H<lo+tp ~ 60 



(?P + ^ P )(?fe + ^fe)-A 2 



-%o-£p + 60 



+ [1 - n(£ p ) - n(&)] 



-%o +C V + 60 



46p6fe 

(6> - Ep)(£k - E k ) ~ A 2 



<%o - 6> - 60 



(£p-£p)(a + £fc) + A 2 



(g p + £ p )fe-ff fc )+A 2 



(Al) 



The axial longitudinal differs only in the sign of the anomalous A 2 term, a consequence of coherence effects being 
different in different channels, and is given by 



Im Il£ {q Q ,q) = -2ir 



d 3 p d 3 k 3 



(2tt) 3 (2tt) 3 



S 3 (p - k - q)(l + p ■ k) 



Hio + tp - 60 



(Cp + E P )(Z k + E k ) + A 2 



-%o-£p + 60 



+ [1 - n(&) - n(60 



(6j ~ E p )(£ k - Ek) + A 2 
8(qo -£ P - 60 



(^p-^p)(6+-Bfc)-A 2 



(C P + Sp)(a-^)-A 2 



46,6c 



(A2) 



From these expressions, one may obtain Iljj , II^' 4 , and H^{ A as specified in Eq. ( p"l| ) 
The transverse response functions have similar forms. The vector is 



Im ir^, (g , g) = -2tt 



d 3 p d 3 k 3 



(2tt) 3 (2tt) 3 



5 3 {p -k- q)(l+p- k - 2p 2 k 2 ) 



Hqo + £ P - 60 



(£p + E p ){£k + Ek) + A 2 



-%o ~ 6> + 60 



+ [1 - n(&) - n(&)] 



-$(qo + 6, + 6c) 



46,6= 

(Cp-^p)(6c-^) + A 2 
46p6c 

%o - 6> - 6c) 



(£p-£p)(6+£fc)-A 2 



46,6c 



(Z P + Ep)(£ k - E k ) - A 2 

46p6o 



(A3) 



and the axial, 
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Im Ilf 2 (q ,q) = -2tt 



d 3 p d 3 k 3 



(2tt) 3 (2tt) 3 



5 3 {p - k - g)(l + p ■ fc - 2p 2 fc 2 ) 



X S K4) - «(6)] 



%o + 6 - 6) 



(4 + £p)(6 + £fc) + A 2 



-%o-6 + 6) 



+ [1 - n(4) - «(&)] 



-5(<?o + 4 + 6) 



466 

(Cp-^p)(6-^) + A 2 
446 

5(go - 4 - 6) 



(ep-^ P )(& + ^)-A 2 



446 



(Cp + gp)fo-^)-A 2 
446 



(A4) 



Finally, there is a small but finite response in the mixed vector-axial channel. Since the neutrino tensor element 
L 23 is imaginary, we take the real part of the polarization. Antisymmetric in spin, it is 



Ren^ A (« R) ,g) = -27r 



d 3 p d 3 k „„ . 

' ■S 3 (p-k-q)(p 1 -k 1 ) 



(2tt) 3 (27t) 2 



x [n(&) -«(&)] 



5(<?o + 4 - 6) 



(4 + £p)(6 + #fc) + A 2 
446 



-5(9o — 4p + €fc) 



+ [1 - n(4) - 



-5(90 + ^ + 6) 



(4-£ P )(4-£fc)+A 2 

5(go - 6 - 6) 



(4-^ P )(6 + ^)-A 2 



446 



(4 + .E p )(4-£; fc )-A 2 
446 



(A5) 
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